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BIFURCATION OF NONLINEAR EIGENVALUES IN PROBLEMS 
WITH ANTILINEAR SYMMETRY 

TOMAS DOHNAL AND PETR SIEGL 


Abstract. Many physical systems can be described by nonlinear eigenvalues 
and bifurcation problems with a linear part that is non-selfadjoint e.g. due 
to the presence of loss and gain. The balance of these effects is reflected in 
an antilinear symmetry, like e.g. the “PT-symmetry, of the problem. Under 
this condition we show that the nonlinear eigenvalues bifurcating from real 
linear eigenvalues remain real and the corresponding nonlinear eigenfunctions 
remain symmetric. The abstract results are applied in a number of physical 
models of Bose-Einstein condensation, nonlinear optics and superconductivity, 
and further numerical analysis is performed. 


1. Introduction 

We consider the nonlinear eigenvalue problem 

Aip - ef(ip) = nxp, (1) 

and analyze the bifurcation in 6 from a simple eigenvalue fi o at £ = 0 in a suitable 
Hilbert space for a rather general class of densely defined, closed (possibly non- 
selfadjoint) operators A and locally Lipschitz continuous nonlinearities /, cf. As¬ 
sumption [(!)] below for details. For a homogeneous nonlinearity /, we also consider 
the additional condition H^H = 1. The main contribution of our paper is to the 
problem of bifurcation from real eigenvalues under an antilinear symmetry of A and 
/. We show that under this condition the nonlinear eigenvalue remains real and 
the eigenfunction remains symmetric. This confirms a number of existing numerical 
computations of specific examples of such a bifurcation problem with an antilinear 
symmetry, see the references below. Besides presenting the abstract bifurcation 
results, we explain in detail how these apply to physically relevant examples by 
checking the assumptions and giving concrete choices of the working space. 

The question of real nonlinear eigenvalues in non-selfadjoint problems with sym¬ 
metries has gained on physical relevance in the recent years due to the intensive 
research on nonlinear systems under the parity and time-reversal (VT) symmetry 
mainly in Bose-Einstein condensates (BECs) [32[ [29] , nonlinear optics |23], see also 
[515 f° r an experimental breakthrough, or superconductivity [48]. In these specific 
physical problems, the presence of real nonlinear eigenvalues typically means the 
existence of stationary solutions of the form with fi £ M also if the system 

is subject to balanced gain and loss (modeled by a non-selfadjoint linear part). 
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The interest in antilinear symmetries was initiated by an observation in [4] 
where Schrodinger operators with VT- symmetric complex potentials in the con¬ 
text of quantum-mechanics-like linear problems were numerically shown to have 
real eigenvalues in a certain parameter region. 

In the context of BECs, where the (nonlinear) Gross-Pitaevskii equation mod¬ 
els the dynamics of the condensate, a complex potential describes the injection 
and removal of particles and a balance of these two processes is reflected in the 
VT- symmetry of the system. Numerical results on the bifurcation of nonlinear 
eigenvalues in particular in one dimensional models can be found e.g. in urnm- 

In optics under the paraxial approximation, the system can be modeled by the 
nonlinear Schrodinger equation (NLS) with a potential corresponding to the refrac¬ 
tive index, which is complex if the amplification and damping of the light wave are 
present, a balance is again reflected in the VT- symmetry. Numerical and formal 
results on the NLS for the bifurcation from linear eigenvalues under “PT-symmetry 
include one dimensional 011571160] or two dimensional 59 . 

Superconducting wires driven with electric currents represent another example of 
a physical application of a PT-symmetric nonlinear eigenvalue problem, cf. [481149] . 
The non-selfadjointness appears due to the dependence of the electric potential on 
the external current. 

As explained in detail in Section [5j our results cover all of the above physical 
models as they are particular cases of 0 with a linear operator A , a nonlinearity 
/ and an antilinear symmetry C compliant with the assumptions of our analysis. 

Our approach to the bifurcation problem 0 is based on the Lyapunov-Schmidt 
reduction and a fixed point iteration. We decompose the Hilbert space to the one 
dimensional ker(A — /io) and its complement using the spectral projection corre¬ 
sponding to the eigenvalue /io and for e > 0 we seek solutions (/i, ?/>) near (/icbV’o), 
where i/jq is the eigenfunction of A corresponding to /io- On the complement of 
ker (A — /iq) the operator A — /i 0 is invertible and a fixed point iteration can be used 
to obtain a small correction of the eigenfunction for 5 small enough. The scalar 
equation on ker (A — /io) is solved likewise by a fixed point iteration and it produces 
a small correction of /io- In this way we obtain an expansion of (i and up to 
second order in e. 

The problem of bifurcation of nonlinear eigenvalues is of course classical and 
has been solved, e.g. in |j] for simple eigenvalues in real Banach spaces and in [28] 
for possibly complex Banach spaces (as relevant in our problem) and for eigenval¬ 
ues of odd algebraic multiplicity or for geometrically simple eigenvalues, see [28l 
Thm.I.3.2]. We choose to prove our results in a Hilbert space and independently of 
[28] in order to provide an explicit expansion of fi and and because the Lyapunov- 
Schmidt reduction and the fixed point equations are used also in the second part 
on the preservation of the realness of /i under an antilinear symmetry condition. 
Moreover, we avoid the technical condition ||/(V ; )|| = ^(IIV 7 II 2 ) as -A 0 of [28] . 

For the problem under the assumption of an antilinear symmetry C of A and /, 
cf. Assumption | (II) [ we check that the fixed point iteration preserves the symmetry 
of the iterates for and the realness of the iterates for fi. As a result, if /io E M 
and if V’o has the antilinear symmetry, then the nonlinear eigenpair (/i, satisfies 
these conditions for e small enough too. 

To our knowledge the only existing mathematically rigorous papers on similar 
bifurcation problems under an antilinear symmetry are [31l and [49] . In [31 the 
concrete example of the discrete NLS with the PT-symmetry is considered. The 
proof is based on the Lyapunov-Schmidt reduction and the implicit function the¬ 
orem. In [49] a one dimensional VT- symmetric nonlinear parabolic problem for 
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superconducting wires is studied using the center manifold analysis. As a special 
case, stationary localized solutions are found. 

The structure of the paper is as follows. Section [2] presents the assumptions on 
the operator A and the nonlinearity / in 0 needed for the general bifurcation 
problem and provides a number of examples of A and / satisfying these conditions, 
whereby we concentrate mainly on Schrodinger operators A but discuss also a first 
order Dirac type operator. Section [2] also explains the choice of our function space in 
which the fixed point iteration is carried out. In Section [3] we prove the bifurcation 
result and the expansion of the eigenvalue and the eigenfunction. For homogeneous 
nonlinearities, we rescale 5 and ip such that a solution with ||^|| = 1 is found. The 
problem under symmetry assumptions is discussed in Section [4j Both antilinear 
and linear symmetries are discussed, where the former one is shown to lead to the 
preservation of the realness of fi. Section [5] explains applications of our results 
to concrete physical problems from literature. Finally, in Section [6] we present 
numerical computations of nonlinear eigenvalues of 0 with A — —A + V and 
f(ip) = \ip\ 2 ip in L 2 (M 2 ) and with VT -symmetric as well as partially VT -symmetric 
potentials V. The effects of a linear symmetry are also observed. 


2. Basic assumption and examples of operators and nonlinearities 

The following basic assumption comprises a condition on a compatibility of the 
linear part A with the nonlinearity / and a spectral condition on A. 

Assumption (I). Let A be a densely defined, closed operator with a non-empty 
resolvent set in a Hilbert space (%,(•,•)) with the induced norm || • ||, let / be a 
mapping in H and let (y, || • \\y) be a Banach space. Suppose that the following 
conditions are satisfied: 

(a) y is a subspace of H, for some n E N is Dom(A n ) C y C Dom(A n_1 ), and 
there are Aq, &2 >0 such that, for all (p E Dom(A n ), 

n —1 n 

\\4>\\n-l ■= < k 2 MVll =: M0||n, (2) 

k =0 k =0 

(b) fi o E C is an isolated simple (i.e. with the algebraic multiplicity one) eigen¬ 
value of A. Moreover, suppose that the normalizations of ip o E Dom(A), 

E Dom(A*), 

Aip o = /ioV’o, A*ipQ = jfiipS, 


i.e. of the eigenvectors of A , A * corresponding to fi o, /io, respectively, are chosen 
such that 

boii = i, (V’o.V’o) = r ( 3 ) 

(c) the mapping / : y —> Dom(A n_1 ) is Lipschitz in a neighborhood of the eigen- 
vector ipo, more precisely: there exist vl > 0 and L > 0 such that, for all 

(t>,4> e{rjey : h-tpoh < r L }, 

\\f(4>) - fWWn-l < hL\\(f> - 1p\\y. (4) 


Remark 2.1 (Remarks on Assumption |(I)| ). The space y is our working space in 
which we perform fixed point iterations. A natural choice for y is (Dom(A), || • 
||i), i.e. the domain of A equipped with its graph norm. Nonetheless, it may be 
convenient to work also with a different y, e.g. with the form-domain and the norm 
induced by the quadratic form of A since these can be much better accessible than 
(Dom(A), || • ||i) itself, cf. Section [2.1| for examples. Obviously, if the Lipschitz 
continuity 0 is established with || • ||y, it holds also with || • || n . A motivation 
for considering n > 1 is given in Examples |2.4| and 2.8, see also Remark 2.9 The 
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condition p(A) 7 ^ 0 guarantees that also Dom(A n ), n > 1 , is dense in 'H, therefore 
also y is dense in T-L. 

Recall that if po is a simple isolated eigenvalue of A , then JIq is a simple isolated 
eigenvalue of A*, cf. [30, Chap.III. 6 .5-6]; moreover, it can be easily verified that 
the normalization © can be achieved. In detail, (^ 0 ,^ 0 ) = 0 implies that V’o £ 
Ker(A —/i 0 ) ClRan(A —/i 0 ) since Ker(A* — po) ± = Ran(A — /i 0 ). The eigenvalue po 
is simple, so Ker (A — po) 2 = Ker(A — po) in particular. From (A — = V’o for 

some (j) G Dom(A) we get (A — po) 2 (t> = 0, thus (j) = expo and hence ^0 = 0, which 
is a contradiction. 

The spectral (Riesz) projection Pq on Ker(A — /io), defined as a contour integral 
for sufficiently small S > 0 , cf. |30j Chap.III. 6 ], and the complementary projection 
Qo := I — Pq can be written, using ([3|, as 

p o~-P [ (A-z)- 1 dz = (;r o )i’o, Qo = I-{; r 0 )%- (5) 

2m JdBsM 

We analyze several groups of operators and nonlinearities below and show that 
they satisfy Assumption | (I) | The selection is inspired by various physical models 
from literature, cf. Section [5j where we apply our results to problems possessing 
typically additional symmetries, cf. Section [4j 

2.1. Schrodinger operators and the space y in Ass. 

operators are naturally associated with the following spaces 

(A II • b) = (U s Q m, || . b & ) := ((H‘(Sl) n Dom(Q)), || • ||*. + ||Q • ||), ( 6 ) 

where Q is a domain in R d , s > 0, Q £ L 2 oc {Q) and Dom(Q) := {ip £ L 2 {Q) : 
Qip £ L 2 (Q)}. It is not difficult to verify that this y is a Banach space. 

Example 2.2 (Schrodinger operators with complex potentials, n = 1). Let 'H = 
L 2 (R d ), V e lE 1 1 0 ’ c 00 (M d ) and V 2 £ L 2 oc (R d ). Let Vi and V 2 satisfy further 

(i) Re Vi >0 and |Wi| < C\ \ V\ \ + C 2 with some C \, C 2 > 0 , 

(ii) there exist a £ [ 0 , 1 ), /? > 0 such that, for all ip £ / Hy 1 (R d ), 

ITaV’H < a(||AV>|| + IIUV’II) + civ’ll* 


( 1 ) 0 . Schrodinger 


Then the operator 

A:=-A + W + R 2 , 


Dom(A) :=U 2 Vl ( 


*) 


and the space y = 'Hy 1 (R d ) satisfy Assumption (I) (jaj with n = 1. The main step 
needed to justify © is the fact that, for all ip £ Dom(A), 


cibllb 


<||A ^|| 2 + b|| 2 <c 2 b||b , 

V 1 V 1 


where ci, c 2 > 0 are independent of ip, cf. for instance [6 j for details; the standard 
proof for V\(x) = x 2 and d = 1 can be found e.g. in |5j Ex.7.2.4]. 


Example 2.3 (Schrodinger operators with singular potentials, n = 1). Let H = 
L 2 ((—r,r)) with r G (0, +oc], V\ G L 1 1 oc ((—r,r)) and v 2 be a sesquilinear form. Let 
Vi and V 2 further satisfy 

(i) ReVi >0 and |ImVi| < tan^ReVi with 0 G [0,7r/2), 

(ii) T~L^/ ReVl ((~A r )) C Dom(r 2 ) and there exist a G [0,1), /3 > 0 such that, for 

all ^ e 7 b'R^d(- r ’ r ))’ 


MV>]| < “(IIV ’ , || 2 + llx/Re Vibll 2 ) + / 5 |lbl 2 - 


( 7 ) 
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Then the m-sectorial operator A associated (via the first representation theorem 
[3UI Thm.VI.2.1]) with the closed sectorial form 

a[ip]:= W\\ 2 + ( Vi(x)\il>(x)\ 2 dx + v 2 [ip], Dom(a) := r)) (8) 

J —r 

and the space y = T~t]y ReVl ((~ r ^ r )) satisfy Assumption I (I) I ([a]) with n = 1. To 
show recall that, for all 'ip G Dom(A) C Dom(a), |a[^]| = \(Ai/j,i/;)\ < (||A0|| 2 + 
||^|| 2 )/2 and the norms || • ||^i_ and W\a[ • ] | + cjj • || 2 with a sufficiently large 

y/Re V 1 

c > 0 are equivalent, cf. [3D1 Chap. VI.]. 

The domain of a and the space y can be selected also, for instance, as 

: = H o((~r,r)) C Dom(yiUU) 

i.e. Dirichlet boundary conditions are imposed at ±r, and the analogues of all 
claims above remain true. 


Example 2.4 (Schrodinger operators with bounded and regular potentials, n > 1). 
We present an example of Schrodinger operators satisfying Assumption | (I) | (|a|) for 
n > 1. The motivation for n > 1 comes from the condition s > d/2 for H s (K d ) 
needed for polynomial nonlinearities in Example |2.8| below, guaranteeing that the 
polynomial nonlinearities satisfy Assumption |(I)|(c 
Let H = L 2 (R d ) and let V G W^(m-i),oo ( 
complex potential. Then the operator 


for some m G N be a possibly 


A:=-A + E, Dom(A) := H 2 (R d ), (9) 

and the space and y = %o n (M d ) = H 2n (R d ) satisfy Assumption |(I)[ (|a|) with any 
n <ra. The claim holds since Dom(A n ) = iL 2n (M d ) and, for all ^ G LP™(M d ), 

where ci, C2 > 0. The second inequality above follows from V G W 2 ^ m_1 ^ 00 (M d ) 
since A n x/> consists of terms TA 21 !/- 7 ' 1 A l2 V^ 2 ... A 2r W- 7n '0 with ik,jk G {0,1} for all 
k = 1,..., n and such that the highest derivative acting on V is of order 2(n — 1). 
Each term can be thus estimated by c||^||#2n. Remaining estimates follow from 
the equivalence of norms on H 2n (R d ), see e.g. [12j Lem.3.7.2]. 


Example 2.5 (Schrodinger operators with quasi-periodic boundary conditions). 
All examples |2.2[|2.3[ and |2.4| can be combined with quasi-periodic boundary condi¬ 
tions on a bounded domain Q C Taking, for simplicity Q = (—r, r) d , r G (0, oo), 
the quasi-periodicity vectors are k G (—7r,7r] d and we can choose 

V = :={iff € n 2 Vl (fl) : il>(x + 2 rej) = e ik ^(x),V^(x + 2rej) = e ik ^V^(x) 

for all j = 1, —, d and all x G dVt s.t. x + 2 rej G 

in Example 1 2.2 [ 

y = r(( _r ’ r )) := U e ' U VRSVl((- r ’ r )) : W) = e ik ip{-r)} 

with k G (—7r,7r] in Example |2.3| (d = 1), and 

y =Hl n ' k (Q) := G H 2n (n) : D m t/>(x + 2 rej) = e lk ^D m ^[x) for all j = 1, ..., d, 

all m G N d s.t. 0 < \m\ < 2n — 1, and for all x G dCl s.t. x + 2 rej G SS2} 

in Example |2.4[ H ere ej is the j-th Euclidean vector and D m = d ... d™ d d . 
Assumption [(I) | (|a|) is still satisfied with these choices. 
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Example 2.6 (Discrete Schrodinger operator). The difference operator on % = 
Dom(A) = C 2N (N G N) defined by 


(Alp) n 


V’n+l + VVi-1 + i7(-l ) n i>n 
< V’n+l - iiV’n 
k -V>n-1 +b'V’n 


for 2 < n < 2N — 1, 
for n = 1 , 
for n = 2 AT, 


( 10 ) 


appears in the modeling of one dimensional optical lattices [31]. Choosing y = %, 
assumption |(I)| (|a|) is obviously satisfied in this finite dimensional case (e.g. with 
the Euclidean norm). 


Example 2.7 (First order Dirac type operator). An example of a physically inter¬ 
esting operator other than a Schrodinger one is 


(-id x -V(x) -k,(x) \ 

V -k{x) i d x -V(x)J 


( 11 ) 


with V, n G L°°(M). This operator occurs in a model for optical waves in fibers with 
a Bragg grating and a localized defect [22] . Assumption (I) 
the natural choice T-L = L 2 

H\\y '= HA\m + \\^2Wh 1 • 


x L 2 


,y = Dom(A) =H 1 


is satisfied under 
x i^ 1 (M), where 


2 . 2 . The spectral condition (Ass. Eq}@)- To satisfy Assumption |(I)| (|b|), a 
detailed spectral analysis of a given linear operator A must be performed. Here 
we recall some perturbation results on isolated eigenvalues that can be often used 
to justify the presence of simple isolated eigenvalues for more complicated, typi¬ 
cally non-self-adjoint, differential operators with complex coefficients. In Section 
[4] Remark |4.3[ we further explain the stability of realness of simple eigenvalues if 
A possesses a certain symmetry. Finally, we recall here basic spectral properties of 
Schrodinger operators, particularly the results on the essential spectrum and the 
simplicity of the ground state. 


2 . 2 . 1 . Holomorphic families of operators. Standard results on the spectrum of a 
holomorphic family of operators A(y), 7 G C, cf [130, Chap.VII], yield that if the 
spectrum of A( 0 ) is separated into two parts, then this remains true also for A(y) 
with |y| sufficiently small, cf [30] Thm.VII.1.7]. Moreover, isolated eigenvalues de¬ 
pend analytically on 7 and their multiplicities are preserved, cf. [30j Thm.VII.1.8]. 
Criteria for the holomorphicity of an operator family can be found in !H Chap.VII]. 
A sufficient condition for the operators or quadratic forms (and hence for the op¬ 
erators associated with a( 7 ) by the first representation theorem [301 Thm.VI.2.1]) 
of the type 

A(y) = A 0 + 7 B, < 2 ( 7 ) = a 0 + 76 , 7 G C, 

where A 0 is a densely defined closable operator and ao is a densely defined clos- 
able sectorial form, is the relative boundedness of H, b with respect to Ao, Reao, 
respectively, i.e. 

A( 7 ) : Dom(Ao) C Dom(H), ||£^|| < a||A 0 ^|| +/3IHI, ^ F Dom(A 0 ), 
a( 7 ) : Dom(a 0 ) C Dom(b), I^MI 2 < aRea 0 [^] + /^I^H 2 , ijj G Dom(a 0 ), 
with some a, /3 > 0, cf. mi Thm.vii. 2 . 6 , Thm.vii. 4 . 8 ]. 


2 . 2 . 2 . Spectra of Schrodinger operators. We consider Schrodinger operators in the 
setting of Examples |2.2f|2.4| many of the following spectral properties are valid in 
a much greater generality, cf. ED] for instance. 

and let lim^i^Qo |Vi(x)| = 00 . Then 


Let A be the operator from Example 


2.2 


the resolvent of A is compact, hence cr(A) = croi SC (A). Moreover, the selfadjoint 
operator A 0 := —A + Q, with Q G L 2 oc (R d ), Q > 0 and lim\ x ^ 00 Q(x) = + 00 , 
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with V = V 1 + V 2 G L° 


has the simple ground state, i.e. the lowest eigenvalue, cf [471 Thm.XIII.47]. For 
d — 1, a Wronskian argument can be used to conclude the simplicity of all eigen¬ 
values for single-well potentials like Aq = —d 2 + |x|^, p > 0. For the singular 
Schrodinger operators from Example 2.3, the resolvent is compact if r < oc (also 
if Dirichlet or quasi-periodic boundary conditions are considered) or if r = oc and 
lim | a .|_ K50 |Vi (x) | = 00 . The simplicity of the ground state in the selfadjoint case 
can be in some situations concluded from [471 Thm.XIII.48]. 

Let A be the Schrodinger operator from Example 1 
and lim\ x \^ cx (Vi + V 2 )(x) = 0, then cr ess (A ) = [0, + 00 ), cf [20, Cor.X.4. 2 , Ex.X.4.3] 
and e.g. 03 Ex.XIII.4. 6 ]. (Note that there are several different definitions of es¬ 
sential spectrum for non-selfadjoint operators, cf. [20, Chap.IX], nevertheless, all 
coincide for this special situation.) Discrete eigenvalues may appear outside essen¬ 
tial spectrum. Particularly in d = 1 , 2 , the selfadjoint operator Ao := —A+sQ with 
f Q dx < 0 and Q decaying sufficiently fast, cf [53j[33i for precise assumptions on Q, 
possesses a unique negative simple eigenvalue for all sufficiently small 6 > 0 ; some 
non-self-adjoint extensions can be found in [45, 46]. For the singular Schrodinger 
operators from Example 2.3 with r = 00 , a ess (A) = [0, +oc) if lim^^co V\{pc) = 0 


and V 2 = 0 or v 2 are forms corresponding to S potentials discussed in Section [522] 
Let A be the Schrodinger operator from Example | 2 . 2 | with V = V\ + V 2 G 


r, r) d ) an d the M-quasi-periodic boundary conditions, k G (— 7 r, 7 r] , so, as 
in Example 


2.5 


Dom(A) = Ho’ fc ((— r, r) d ). These operators arise naturally from 
a periodic problem in L 2 (M d ), where V is extended periodically onto cf 071 


Chap.XIII. 16] and the concepts of the “band structure” and of the Bloch eigenvalue 
problem. As mentioned above, the resolvent of A is compact. Moreover, if d = 1 
and V is real, all eigenvalues of A are simple if k £ { 0 , 7 r}, cf [47, Thm.XIII.89]. 


2.2.3. Spectrum of the discrete Schrodinger operator (10). It is a straightforward 


calculation, see m, to show that A in (flQ[) has the 2 N eigenvalues 


(■ 


± 4 cos 


2 ( nj \ 


1 


1/2 


l<j<N , 


\l+2N J 7 

which are simple and real for 7 G ^—2 cos ( 37 ^ 7 ) 5 2 cos ( jqpf/v )) • 

2.2.4. Spectrum of the Dirac type operator jnj. if k(x) -a /Coo > 0 and V(x) — » 0 
as \x\ oo, then the essential spectrum of A in ( [IT] ) is (—00, —k qq] U [n oo, 00). To 
see this, notice the perturbation result uni Prop. 6.6] and the unitary equivalence of 
A in ( 11 ) and H in [TO] for \x\ -A oc. For real n and V the operator is self-adjoint. 
Bounds on eigenvalues of non-self-adjoint perturbations outside of the essential 
spectrum are proved in m- Special choices of real n and V with simple eigenvalues 
in (—/Coo, ^ 00 ) have been found in (22] Sec.4.B.] using a connection to the spectral 
problem of the inverse scattering theory for the nonlinear Schrodinger equation. 
For example, for k(x) = (jliq + k 2 tanh 2 (/cx )) 1 / 2 and V(x) = k~ 2 (x)sech 2 {kx) 

with //o G (—/c 0 
Mo- 


, /Coo) and k G M \ {0}, the operator A has a simple eigenvalue at 


2.3. Nonlinearities (Ass. [(I)j (|c|)). We present several nonlinearities / satisfying 
Assumption HI©- The considered spaces y are those arising for Schrodinger 
operators in Section |2.1| In fact, we show that the Lipschitz continuity © holds 
locally for all 770 G T, thus also for the eigenvector 7 /V 
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Example 2.8 (Polynomial nonlinearity). Let % = I/ 2 (fl) and y be as in |6| with 
s > d/2, sgN, and £2 = R d . Let iVeN and 


N 


/poiWO := a pq W with a pq £ C^(R d , C) for all p,q , 
p,q =o 


( 12 ) 


where C§ is the space of functions with continuous and bounded derivatives up 
to order s. Without any loss of generality, we set aoo = uoi = aio = 0. A 
classical example is the cubic nonlinearity / c (0) := |'0| 2 '0, ^.e. a 2 i = 1, = 0 

otherwise. We show below that / po i in (12) satisfies Assumption E30 with any 
Vo = n s Q (R d ) and tl > 0. 

First note that, for s > d/2, the norm || • || h s satisfies the so-called algebra 
property: there exists C a > 0 such that, for all 0, 0 G H s 


WHWhs <^8,11011^11^11^, (13) 

cf. [U Thm.4.39] or [T7l Lem.4.2]. Moreover, the Sobolev embedding of H S {£1) 
in holds, cf. [1, Thm.4.12], i.e. there exists C c > 0 such that, for all 0 £ 

H s (R d ), 

\\4>\\l°° < C e \\<f)\\H‘- 

Thus, for all (f),ip £ H s Q (R d ), 

Ill’ll H% = WHW h- + WQHW 


(14) 


< Ca,\\(t>\\H°\\4’\\H° + VMl-Ml-IMW iiedi 

< maxjca, UW-H^Mln^y 

hence the norm || • \\^s satisfies the algebra property as well. 

Clearly, it suffices to check @ for a single term From 

- vr = ^ p )(^ q +r) + \(<t> p +v)$ q - r), 


(15) 


m— 1 


4> m _ = (4>- if) 


-1 -k 


(16) 


k =0 


and using (15), we obtain 


wr - vru h < ci (n^ - vu h w + ru h + w+r\w h w - ru h ) 


< c 2 


k =0 
q -1 


m p n h + mc & ) Eii^ii^iiS 


k =0 


-1 -k 
Q 


II U; 


Q’ 


where > 0 are independent of 0 and 0. Finally, for all £ G {v G Hg( 

< r L } with any % 6 ?fg(K d ) and r L > 0, we have ||£||«^ < ||r? 0 ||«^+r-L, 

thus 

wr - rr\\u h < c v 0 , rL , P , q u - nu h 


if 4>, ip £ {q £ n s Q 
f,£W Q 


— t)()\[h q < r L }• Also note that for a £ C^(R d ) and 


\\ak\\n‘ Q < IMIc s II^IN s + IMIc 0 IIQV’II < IMIcHIV’lkj,- 
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As a result 


|J/pol(0) - /polWOllw, < 



H - i>\\w c 


(17) 


which implies validity of 0 in Assumption | (I) | (|c|). 

Note that this Lipschitz continuity can be directly extended to vector valued 
polynomial nonlinearities ffy) = (/i(V0> • • • ? /m(V0) T with m G N, 7 / G (%g( 

and with fj(ifi) = E 


AT 


,09 


., Prn ,q rn =0 Uj PlQl---Prnqrn' l l J l 1 ' l Pl * * * Vfr V^? 
must be replaced, e.g., by := Ejli ll^illw 


Qrri 


case the Hg-norm in 


o 


In the vector 


Remark 2.9 (Motivation for n > 1 in Ass. [(I)| ([a|)). Examples |2.4| and |2.8| constitute 
the primary motivation for n > 1 in the choice of a spa ce y with Dom(A n ) C^C 
Dom(A n_1 ). The condition N 3 s > d/2 in Exam ple |2.8| implies that for d > 4 
we need s > 3. The natural H 2 -space of Example |2.2| is, therefore, not sui table 
as our working space. On the other hand, the space H 2n of Example |2.4| with 
n > d/4 is sufficiently small. The reason for choosing Dom(A n ) C y C Dom(A n_1 ) 
is the use of a fixed point argument in Theorem in particular to guarantee 
G(B r2£ 2 ) c QoE cf. the reasoning between (26) and 
In summary, if we choose 


d s 

- < - < n < m with s,n,m G N, 
4 2 _ _ 


(18) 


then A in ( j9| w ith V G W 2 ( m 1 ^°°(M d ), / = / po 1 , and y = id 2 n (M/) satisfy 




2.5 


with the 


Assumption (I) 10 and 0- Of course, also y = H 2n,k (Q) of Example [ 
same A and / is admissible. For d = 1,2 inequality Jl 8 |) holds also with 8 = 2, n = 1, 
such that y ■ 


of Example 2.2 or the space of Example 


2.5 


can be 


used with / = f po \. For d = 1 we can use eve n s = 1, n m 1 and, hence, the spaces 
((—r, r)),H ^ eVi ((—r, r)) of Example 2.3 and 4^ ^ ^ ((—r, r)) of Example 


Re V\ v 


12.51 are admissible. 


VReUi V 


Example 2.10 (The monopolar and dipolar interaction from [38l Hi] ). Let T-L = 
L 2 ( 0 ) and 3 ^ be as in 0 with 8 = 2 and Vt = M 3 . We investigate nonlinearities 
having formally the form 

(f L {^))(x) \=^{x) [ K L (y)\^{x-y)\ 2 dy, ^G{m,d}, 

J R3 

with kernels K m (x) = and i^d(^) = jyp ^1 — 3^ with fixed a G M 3 , 

||a || R 3 = 1 . Notice that choosing id(x) = d 3 (x), we can recover also the so called 
contact interaction, which, however, coincides with the cubic nonlinearity from 
Example |2.8| Here we focus on the nonlinearities / m and /d, the so-called monopolar 
and dipolar interaction, respectively, see Section pu 2 | for more details. 

First we show that, for any 77 G H 2 (M . 3 ) fi L 1 (M 3 ), 


/ 

J R ! 


rj(x - y) d y 


\v\ 


< C\\t)\\ H 2 + IMIli 


L °o 


with C independent of 77 . Indeed, the Sobolev embedding (14), applied in the last 
step, yields 


/ 

J R 3 


\v(x - y )I d y 


-L 


\v( x ~y)\ dy 


+ 


B i(0) 


/ 

d B\ (0) 


/ 

iR 3 \Bi(0) 


\t](x -y)\dy 


dy 

\y\ 


+ IMIli < GMI # 2 + IIjjIIli. 
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Using this, the algebra property of || • ||# 2 , cf. (13), the special case of the first 
formula in © and Cauchy-Schwartz inequality, we obtain 

\4>{x - y)? d y 


\\fm(4>) - fn 


<u~- 


+ 


I 

JR 3 


/ 

J R 3 


\y\ 

(x - y) I 2 - I ip{x - y )| 2 | Ay 


\y\ 


< Ci\\ 4>\\ 2 H 2\\4> --ip\\ H 2 Q + c 2 Ill’ll \\4> — V , IIh*IP + V'IIh 2 


\\h 


<c 3 m\h + 

thus © is satisfied with n — 1 for any 770 £ %q(M 3 ) and El > 0. In particular, 
Assumption | (I)| (|c|) holds. 

The nonlinearity /<* is more complicated and it is even not immediately clear 
why it is well-defined. Nonetheless, the appropriate framework of singular integrals, 
particularly [54, Thm.II.3], yields that, for any 77 E L 2( 


I 

Jr 3 


Kd(y)v( x ~y) d y 


< CdM, 


with Cd independent of r\. Using the Sobolev embedding ([14]), the algebra property 
(13) of iL 2 -norm and the special case of the first formula in (16), we obtain 

\\m)-f d m<c 4 m 2 i 


is 2 


+ 


iff 2 ) 


thus Q is satisfied with n = 1 for any rjo E Hq(M 3 ), vl >0, and Assumption 
|(I)|(|c|) holds in particular. 


Example 2.11 (Nonlocal nonlinearity from [49]). Let T-L = L 2 (( — 1,1)) and y = 
Hj’ D ((-l, 1)) = 1,1)). We analyze the non-local nonlinearity 

(/nW)W •= / Im (■ ip(s)ip x (s )) ds. 

Jo 

The validity of 0 is easily checked since /n satisfies 

\\h(4>) - /nWOII 


nX nX 

(cp-ip) Im (<f4> x + ip'il’x) As + + 'ip) / Im (<M> X - 4>il> x ) ds 

Jo Jo 

(||0-^IIl~ J o 


< 2 ( U-^Wl- 
+ \\ 4 > + V’llj 


/ 


Im [4 
Im 


+ V’V’*) ds 

- ds 


For any r) 0 e i?o((-l,l)), r L > 0 and (f>,ip E {rj E : h-yoWm < x l}, 

we have 


f 


Im 


+ i’i’x) ds 


< 


>11+11^1111^x1111 < 


i + HV’llffi) 


< V2(r L + hollffO 


and 


f 


Im 


- x/jip x ) ds 


1 

< - 
~ 2 


' Vh 


< v / 2||0 + V ; l|R rl ll0-' 

< 2a/2 {r L + ||?7o||tti)|| 


I// 1 




I# 1 - 


'Ipx 
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Thus, by the embedding of H 1 in L°°((— 1,1)) as in Example 2.8, we obtain 0 - 
In summary, Assumption | (I) | (|c|) holds with n = 1. 

3. Nonlinear eigenvalue problem 

First, we prove the local existence and uniqueness of the solutions of the nonlinear 
eigenvalue problem 0 under Assumption m Next, we focus on homogeneous 
nonlinearities, for which 0 together with condition \\i/j\\ = 1 can be solved. Finally, 
the influence of possible non-selfadjointness (more precisely non-normality) of A on 
constants appearing in our estimates is discussed. 

3.1. Local existence and uniqueness. 

Theorem 3.1. Let A and f satisfy Assumption El Then every solution of the 
nonlinear eigenvalue problem 

(A - - ef(ip) = 0, (V’,V’o) = 1 (19) 

can be written as 

M = Ho + £V + £ 2 <J, ip = Ipo + £(/) + X, (20) 

with u, a £ C and <j>, x G QoDom(A) and where 

v = -(/(V’o) ,r 0 ), ( 21 ) 

(j) is the unique (in QoDom(A)) solution of 

(A - fi 0 )4> = vtpQ + /(V’o), (22) 

and (a, x) solves the nonlinear system 

0 = £<r+(/(V>)-/(V’o),V’o). ( 23 ) 

Qo(A - tio)QoX = £[(" + ecr)(x + £<t>) + <9o(/(V0 - /(V’o))] =: R{x)- (24) 

Moreover, there exists > 0 such that, for all £ E (— 60 , £ 0 ), the nonlinear 


eigenvalue problem (23)-(24) has a unique small solution, namely with 

ji 


W\ < r i , 


|y < r 2 e 


where r\,r 2 = 0(1) (£ —> 0) satisfy (28), (32) respectively. 


Proof. Without any loss of generality, for a given (pi 0 , V’o) we can write the solution 
(/i,-0) as in ( |20| ). Although at this point the representation is not unique, we 
already know that £f> + x = Qo(£<fi + x) due to the constraint (-0, ^ 0 ) = 1 and the 
normalization (^o^o) = 1- 

First, we apply the projection Po to ( fl9| ) and obtain 

v + sv= -(/W 0 ,^o) = “(/(^0)^0) “ (/WO - /W>o),^S)- 


Defining u := — (/(V’o), V’o)? we get equation (23). 

Second, we apply Qo to ©, resulting in 

Qo(^ - po)QoN> + x) =0^ + e 2 (y)(e(f + x) + eQo/WO 

=£ [(v + 6Cr)(50 + X) + Qo(/W0 ~ /W>o)) + W’O + /W>o)] • 

Next, we define <p to be the unique solution E QoDom(A) of ( [22] ) . This solution 
exists because + /(V’o) £ Ker(A* — po) -1 = span(W)}“ L and because /io is not 
in the essential spectrum cr e 5 (A), it is not in cr e s(A) either, cf. [20j Chap.IX], and 
hence A — pi 0 is Fredholm. The equation above thus becomes problem (24). 


The rest of the proof deals with the existence of a unique solution (cr, y), with x 
small and a bounded, of (23)—(24). Note that Qo(A — /io)Qo is boundedly invertible 
in QpH , hence equation (24) can be rewritten as 

X = (Qo(A- f i 0 )Qor l R(x)=-G(x )• (25) 
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In the first step, for 5 in a small neighborhood of 0, we use the fixed point 
argument to conclude the existence of a solution x of (25) with ||x|| y = 0(s 2 ). We 
search for a fixed point 

x e B r2S 2 := {q e Q 0 y : ||j?|| y < r 2 e 2 } (26) 


with some > 0 independent of 6. Its existence is guaranteed if we can show 

(i) X £ B r2£ 2 => G(x) £ B r 2 £ 2 , 

(ii) there exists p G (0,1) such that ||G(xi) — G(x 2 )\\y < p||xi — X2II3; for all 
Xi? X 2 £ B r2 s 2 - 

Note that ^ 0 , being an eigenfunction, satisfies i/jq G Dom(A m ) for any m G N. 
Thus the right hand side of ( [22] ) lies in Dom(A n_1 ), hence 0 G QoDom^ 72 ) c Q^y. 
Therefore if % G then Qo/(V0 G QoDom(A n_1 ) and R(x) £ QoDom(A n_1 ). 
It is straightforward to check that then G(x) £ QoDom(A n ). The properties of the 
norm || • ||y, cf. Assumption | (I) | (|a|) , yield 

l|G(x)b < |l|G(x)IU < ^C, 0 \\R( X )\\n-i- (27) 

To show the second inequality in ([27]), note that because G(x) ^ QoDom(A n ), 

AG(x) = (A — /io)QoG(x) + Mo G(x) 

= Qo(A — mo)QoG(x) + Mo G(x) = R(x) + Mo G(x)- 


We thus have, for k > 1, 


TG(x)=T- 1 (i?(x)+MoG(x)) = ' 


fe-1 


i=0 


Asaresult, ||A fe G(x)|| < Ej=o iMol* 1 i ||Ti?(x)|| + |Mo| fc ||(Qo(A-^ 0 )Qo) 1 |lll i 7(x)ll 
and the second inequality in ( [27] ) follows, where C ^ 0 > 0 is a constant depending 
on |Mo|,rc and \\(Qo(A - Mo)Qo) _1 ||- 

To ensure ([!]), we take x £ P r e 2 > estimate ||P(x)|| n -i and select suitable r in 
the following. First note that Po and Qo are bounded on Dom(A m ), m G No, 
moreover, since A k Po!p = P^A k ^ and A k Qo'ip = Q^A k ip for all ijj G Dom(A m ) and 


fc = 0,1,..., m, we have 
\\Po\\m := sup 

0 7 ^ p G D o m ( A m ) 


II-PqV’IU 

Ibllm 


sup 

0^</>eDom(A m ) 


Er= 0 iia^vii 

Er= 0 M'bii 


< IToll, 


and similarly ||Qo||m < ||Qo||- Next, 

|l-R(x)||n-i < |e|(k + ^|(||x||„-i + NblU-i) + !]<2o(/(V>) - /bo))IU-i) 

< |e|fei [\v + £a\(r 2 e 2 + |e|b||y) + ||Qo||-P||e<)> + xllx) 

< e 2 k\(\v + e<j|(r 2 |e| + ||</>b) + \\Qo\\L(\\<t>\\y + r 2 \e\)), 
thus we select r 2 such that 


r 2 >k 2 C„ 0 (\is\ + \\Qo\\L)\m y . (28) 

For all sufficiently small 6, we satisfy firstly H^ — ipo\\y < kl W^Wy +^ 2 ^ 2 < Ml, hence 
the Lipschitz property of /, cf. Assumption EH©. can be indeed used. Secondly 
we satisfy condition di ). 

It remains to prove (m. Similarly as above, we obtain (with 7/^2 := ^o+^+Xi^) 
ll^(Xi) -^(X 2 )||n-i = \e\\\(y + scj)(xi ~ X 2 ) + <3o(/(Vh) - /W> 2 ))||n-i 
< |e|fci(|i/ + ecr| + ||Qo||A) ||xi — X 2 1|A 2 * 

Hence, for all sufficiently small 6, condition ([h]) is also satisfied. 
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In summary, there exists So > 0, such that, for all 5, |e:| < £q, we have the 
function x £ ^r 2 e 2 that solves (25); note that then x £ QoDom(A) as well. Note 
also that x and i n particular £q depend on a. However, we consider only \a\ < rq, 
where rq satisfies (32), and an inspection of the estimates above shows that we can 
find £o independent of a (dependent only on rq). 

In order to solve the first equation in (23), we prove first that the solution x is 
continuous in a. More precisely, 


HxOi) - xMb < |qj|i?(xUi)) - R(x(^))\\n-i 

= I^JelllKxUi) - x(ct 2 )) +e((Tix(cri) - <t 2 x(<7 2 )) + e 2 {?\ ~ o 2 )<j> 
+ Qo(f(lpO +£(/) + x( cr l)) — /Oo + s4> + x(CT 2 )))||n-l 

< ^CJ^MlIxUi) -xU 2 )||y +e 2 \(J 1 ~ a 2 \U\\y 

+ y IIUi - o- 2 )(x(o-i) +x(ct 2 )) + (cti +ct 2 )(x(oi) -xU 2 ))||y 
+ IIQoll^UxUi) -x(o- 2 )||x) 

< hCf, 0 |e| ^|e|(I||xUi) +x(o- 2 )||y + NII0||y)ki - a 2 \ 

+ (M + IIQolli + y Wi + c 2 |) ||x( eT i) _ xMb ) ) 


Hence, for \cfi^\ < n, 

HxUi) - xU 2 )b(l - k 2 C^ 0 |e|(M + \e\ri + ||Qo||L)) 

,, x (30) 

< k 2 C IXo \e\ 3 (r 2 \e\ + \\4>\\y)\cn - o 2 1 . 

As the final step, we use the fixed point argument on 

O’ = -I (/(VO - /(V’olV’o) = : S(a), (31) 

where we search for a fixed point in {a G C : |cr| < rq} with a suitable rq selected 
below. Since 


we choose rq such that 


|5(a)|< fc 1 L||A||(bb + r 2 |e|), 


ri >hL ||P 0 ||My> 


(32) 


hence, for sufficiently small |e:|, \a\ < rq implies | S(a)\ < rq. Moreover, using the 
continuity of x m c f • (30), we obtain 


| 5 (<ji) - S(<t 2 )| < ||/(V’Ui)) - /(V’(o' 2 ))|| 

< yrllxUi) - x(o- 2 )b < Ce 2 \a x - <r 2 |, 


hence, for sufficiently small |e|, the fixed point argument yields the sought solution 
of pll). □ 
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3.2. Homogeneous nonlinearity. In the case of a homogeneous nonlinearity like 
e.g. f(if) = | if\ 9-1 V’? solutions with norm one can be generated from the nonlinear 
eigenfunctions of Theorem |3.1| by a scaling. 

Corollary 3.2 (Nonlinear eigenfunction with norm one). Let A and f satisfy As¬ 
sumption [(I)] and suppose that f is a homogeneous nonlinearity, i.e. for all a > 0, it 
satisfies the scaling property f(aip) = a q f(if) with some q E M. Given the nonlinear 
eigenpair (/i(e),^(e)) for |e:| < £q from Theorem \3.1\ the pair 

(MAMA = (Me), HV’OOirVOO) withe = e ||^(e)|| 9_1 

solves (A — — if ('if) = 0 and satisfies ||'0(£ r )|| = 1- The mapping £ t-A £|| / 0(£)|| g,_1 

is injective for \e\ < £\, where £\ < £o is small enough. 

Proof From the scaling property we immediately get 

(a - Me)) nm-'m - = o. 

The injectivity follows from the asymptotic equivalence e\\if(e) H 9-1 = £\\ifo + £(f + 
x\\ q ~ X ^ e||^o|| 9_1 = e for 5 -a 0. □ 


Remark 3.3. Similar results on the bifurcation from simple eigenvalues of Fred¬ 
holm operators appear in the literature. As mentioned in the introduction, a clas¬ 
sical reference is the paper |28j by Ize, where Thru.1.3.2 applies under our assump¬ 
tions and the additional condition ||/(V ; )|| = OdlV’H 2 ) as if -A 0. The theorem of 
Ize treats the bifurcation from the zero solution if = 0 in (A — pi) if = f{if) but 
for a homogeneous / the problem can be rescaled to §• Our Theorem |3 .1 1 avoids 
the technical condition ||/(V ; )|| = Odl^ll 2 ) idf ~^ 0) and provides a more explicit 
expansion of p and if. 


3.3. Constants and non-selfadjointness. Our fixed point argument in Theorem 
|3.1| works for £ sufficiently small and the size of remainders x and c is determined 
by constants r<i and rq, respectively, cf (28), (32). Notice that the size of e is 
restricted at least by 


(W + + ||Qo||i) < 1? 


arising from the contraction condition p < 1 together with ( |27| ) and (29). 

Provided A is normal, i.e. AA* = A*A, the spectral projection P 0 and the 
complementary projection Qo are orthogonal, hence |jPo|| = ||Qo|| = 1, and, in the 
case n = 1, the constant is determined by spectral properties of A since 


C/io — 1 


(1 + \po\)\\(Qo(A - po)Qo) X || = 1 ■ 


1 + I Mo | 


dist(/i 0 , cr(A) \ {/x 0 })’ 


due to the standard relation || (A — z) _1 || = dist(z, a( A)) -1 valid for normal oper¬ 
ators. 

However, in our applications, we usually encounter non-symmetric perturbations 
of self-adjoint operators that result both in non-selfadjointness and non-normality 
of the perturbed operators. Note that we lose both equalities for ||Qo|j and if 
A is not normal, only > is left in general. In particular, the spectral projections as 
well as the complementary projections may behave wildly as the size of the spectral 
parameter is increased even for simple looking one-dimensional Schrodinger opera¬ 
tors with a complex potential and compact resolvent. For instance, considering the 
rotated oscillator —d 2 + ix 2 , cf. [T3l fl4] , for which all eigenvalues A n are explicit, 
X n = e 17r / 4 (2n +1), n E No, the norms of the corresponding spectral projections P n 
grow exponentially, more precisely, 


log IIAII 


lim 


= const ; 







BIFURCATION WITH ANTILINEAR SYMMETRY 


15 


similar behavior is exhibited also by other well-studied (often PT-symmetric and 
with real spectrum) Schrodinger operators, cf. for instance [26, [27] 03], 36. and 
references therein. Notice that the growth of ||P n || implies the growth of ||Q n || 
since ||Q n || > ||P n || — 1. Concerning the size of C Mo , the norm of the resolvent of a 
non-normal operator, i.e. its pseudospectrum, see [55], may be dramatically larger 
than dist(z, a(A))~ 1 . While there is a collection of recent pseudospectral results 
for non-normal differential operators, cf. [13113 [S5J, the estimates on the norm of 
the resolvent of (A — fio) \ QoDom(A) acting in QoH seem not to be available. 

On the other hand, there exists a large collection of perturbation results, partic¬ 
ularly for operators with compact resolvent, cf. for instance the classical ffl CMol 
or more recent [5UES1E1I3S1], guaranteeing that the eigensystem of a perturbed 
selfadjoint (or normal) operator contains a Riesz basis, hence ||P n ||, ||Q n || and 
\\(Q n {A — fi n )Qn)~ 1 \\ are uniformly bounded (||P n || are uniformly bounded already 
if there is only a basis). The Riesz basis property is present for example for opera¬ 
tors in ( [41] ), ( [42] ) , (43) and (50) from Section [5] 


4. The role of symmetries 

We show that under antilinear symmetry assumptions on A and /, cf Assump¬ 
tion [(IT)] below, the nonlinear eigenvalue fi starting from a real eigenvalue /io remains 
real and a certain symmetry of the solution ip is preserved for all 6 E (— £q,£q). The 
symmetry of solutions is preserved also in the case of linear symmetries that are 
studied next. 


4.1. Antilinear symmetries. 

Assumption (II) (Antilinear symmetries of A and /). Let A be a densely defined 
and closed operator in a Hilbert space T-L and let C be an antilinear , isometric 
and involutive operator, i.e. for all 0, 'ip E T-L and A E C, C(\(p + 'ip) = A C(p + C'lp, 
(Ccp,Cip) = {ip, <p) and C 2 = /, such that 

(a) for all rp E Dom(A), 

Cip E Dom(A) and AC'ip = CA'ip , (33) 

(b) for all ^ £ y, where y is the space from Assumption | (I) [ (|c|) , 

cm = f(C VO. (34) 

The operator C is referred to as the antilinear symmetry of A and /. The 
standard example of C is the VT symmetry that is naturally present in various 
physical models as we indicate in the following example and in Section [5] 

Example 4.1 (VT symmetry). Let T-L = T 2 (Q), where Q = or Q = (— r, r) d 
with r E (0,+oo). Define 

(V^)(x) := 'ip(-x), T'lp := 'ip. 


In quantum mechanics, V corresponds to the space reflection (parity) and T is the 
time-reversal. The antilinear VT symmetry is the composition, i.e. C := VT. In 
more dimensional domains, the so called partial ViT symmetries, where 

(Viip)(x i,..., x u ..., x d ) := ^(x !,..., -Xi,..., x d ), 


are sometimes considered, cf. [7 or m • Obviously, also ViT is antilinear. 

Schrodinger operators — A + V with complex potentials V, cf. Examples |2.2| 
and |2.4[ are PT-symmetric if V is PT-symmetric, i.e. [VT,V] = 0, or, in other 
words, the real and imaginary parts of V satisfy (ReE)(— x) = (ReE)(x) and 
(lmV)(—x) = —(lmV)(x). The Schrodinger operators from Example |2.3| posses 
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this symmetry if, for every <p,ip G Dom(a), (</>, P7T/0 
Section [5] for examples of such . 

The polynomial nonlinearity / po i(^) from Example 
|(II)| ([b|) with C = VT if and only if a pq (—x) = a^(x) and the space y can be 
sele cted as in Example | 2 . 2 [ |2.3[ |2.4| or |2.5| For the quasi-periodic case in Example 


24 


V 2 ('PT 4 >,'ip) holds, see 
satisfies Assumption 


245] note that ip is /c-quasi-periodic with a given vector k G (— 7 r, 7 r] d if and only if 
'ip(-x) is fc-quasi-periodic. The monopolar and dipolar interactions / m and /d from 


Example | 2 . 1 Q| and the nonlocal nonlinearity in Example 2.11 are also VT -symmetric 
as it can be easily checked. 

We recall simple facts about spectral properties of C-symmetric operators. 

Lemma 4.2. Let A satisfy Assumption \(ll)\ and let /io be an eigenvalue of A. Then 
(i) ~pT is an eigenvalue of A and if there is a C-symmetric eigenvector ^0 corre¬ 


sponding to fiQ, i.e. 


(A - fi 0 )ipo = 0 , Cipo = V’o, 


(35) 


then fi 0 G M. Moreover, if /jlq is real and simple, then the corresponding 
eigenvector can be chosen C-symmetric. 

(ii) if fa 0 is isolated simple and real, then, in addition, the spectral (Riesz) projec¬ 
tion P 0 of A corresponding to /io commutes with C, i.e. 


PoC = CPn. 


(36) 


Proof. (i) Let ^0 £ Dom(A), ^0 7 ^ 0 satisfy (A — /io)ipo = 0- Clearly, by the 
C-symmetry of A, AC'ip 0 = CA'ip 0 = VPCipo, hence ~pA is an eigenvalue of A. 
Next, let ipo be as in (35). Then 


/ioV’o = ACCipo = CAi/jq = C/ioV’o = Mo^o, 


thus /io G M. 

Finally, let again ^0 satisfy (A 


/io)V ; o = 0 and /io be simple. It follows 
from (33) that ACfi 0 = fioCipo, hence both of ipo d= C^o are eigenvectors of A. 
However, one of them must be 0 since /io is simple. Note that ^0 + C^o is 
automatically C-symmetric and if it is zero, then we take i(^o — C^o)- 
(ii) Since A commutes with C, we have C(A — z)~ x = (A — z)~ 1 C. Hence, for the 
spectral projection Po, cf. we get 

1 r 27r 

CP 0 = —C— / (A - (ji 0 + 

^ JO 

{A - (no + 5e~ i, ' p ))~ 1 5e~ 1(p dipC 


1 f 

2tt J 0 

T f° ( 
2tt y_ 2 4 


(A-zy'dzC = P 0 C. 


□ 


Remark 4.3 (Stability of real simple eigenvalues). For a C-symmetric family A(y), 
cf. Section isolated simple real eigenvalues of A( 0 ) remain isolated simple 

and real for small |y| since eigenvalues are analytic in 7 and always form complex 
conjugated pairs, cf. Lemma 4.2 (ij. 

Notice that many examples in literature as well as in Section [5] can be in fact 
viewed as a holomorphic family with A(0) = A(0)*, thus cr(A(0)) C M. A typical 
behavior of real eigenvalues as 7 is increased, i.e. the non-symmetric part of the 
operator becomes stronger, is a tendency to merge and create a complex conjugated 
pair, see e.g. Figure [4] (a), (b). 
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Finally, we remark that there are also C-symmetric (actually with C = VT) oper¬ 
ators with real spectrum that are not small perturbations of a selfadjoint operator, 
e.g. the celebrated imaginary cubic oscillator —d% + ix 3 acting in P 2 (M), cf. mm- 


The following theorem shows that real simple eigenvalues /io of an operator A 
persist to be real nonlinear eigenvalues n(s) of A — ef for all 5 small enough if 
A and / possess an antilinear symmetry. Moreover, starting with a C-symmetric 
linear eigenfunctions ipo of A (the existence of which is guaranteed by Lemma 4.2), 
the nonlinear eigenfunctions 'ip(s) are also C-symmetric. 

Theorem 4.4. Let A and f satisfy Assumptions |(I)| and \(LV)\ Suppose in addition 
that fi o Gl and choose the corresponding eigenvector ipo as C-symmetric, i.e. Cip o = 
ipo. Then, for all £ G (— £o,£o), the nonlinear eigenpair (/x, -0) from Theorem 3.1 
satisfies fi G M and Cip = ip. 

Proof. Recall that /i = po + eu + s 2 a, ip = ipo + scj) + x and Po is the spectral 
projection of A corresponding to the eigenvalue fi q. 

In the first step, we show that v is real and C(p = <p. The spectral projection 
Po can be written as Po = (-,ipo)ipOi where ip^ as m Theorem 3.1 therefore 
—vipo = P 0 (/(<M), c f d^T] )- Using symmetries ( [34] ) , ( [35] ) and (36), we obtain 

(37) 


C(-^o) = CP 0 (/(^o)) = Po{f{Wo)) = Pomo)) = -W’o, 


thus — vipo = —vipo, hence v G 
tries of A , / and ipo, cf. ( 33 ) 


. Applying C to equation (22), using the symme- 
and ( |35| ), we get 

(A - /i O )C0 = vip 0 + /(V’o)- 

([36]) and because the solution of (22) with Po<p = 0 is 


Because PoCcp = CPo<p by 
unique, we have C(p = (p. 

Let us now work on a. Since a is the solution of the fixed point problem a = S(a) 
with S(a) = -\{f{ip{cr)) - f{ip o)Ao)> where ip (a) = ip 0 + £(p + x(cr) and x solves 

the fixed point equation x = G(x;cr), it remains to show that the coupled fixed 

point problem preserves the realness of a and the C-symmetry of y. 

Given a G M (with \a\ < rq), we prove that 

Cx = X CG(x;cr) = G(x;cr). (38) 

As G(x;cr) = (Qo(A — /io)Qo) _1 P(x; <r), we first show the analogous property for 
R and then the commutation of ( Qo(A — /io)Qo) _1 with C. Since Q o = I — Po, we 
get from ( [36] ) that Qo C = C Qo as well. Also note that for Cx = X the full solution 
ip = ipo + f> + x is C-symmetric. Hence, for 6, a G M and Cx = X> 

C-R(x; o’) = e({v + ecr)(x + e<p) + Q 0 C(f(ip) - /(V’o))) 

= e((r / + ecr)(x + £(p) + Qo(f(Cip) ~ f(Cip 0 ))) = #(x; A 


To prove (38), it remains to show that C(Qo{A — po)Qo) 1 = ( Qo(A — po)Qo) 1 C. 
To this end, take any cp G QoTL, then 

C(Qo(A - /i 0 )Qo) _ V 

= ( Qo(A — po)Qo)~ 1 (Qo(A — po)Qo)C(Qo{A — po)Qo)~ 1( P (39) 
= (Qo(A - p 0 )Qo)~ 1 Cp. 

Property ( [38] ) implies that the fixed point of x = G(x) in P r2 £ 2 > c f • (26), lies in 
P r2£ 2 D {77 G P : Ct? = 77}. 

Finally, we need to show that 
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Once again, because C\ = X implies C0 = ^, we get by a straightforward manipu¬ 
lation analogous to ( |37| ), 

C(-e5(cr)0o) = CP o (/(0) - /(0o)) = P o (/(C0) - /(C0 O )) = -eS(a)^ 0 , 

hence S(a) G M by the same arguments as below (37). □ 


4.2. Linear symmetries. The operator A and the nonlinearity / may possess also 
a linear symmetry 5. Then (for simple eigenvalues) symmetry or antisymmetry of 
the nonlinear eigenfunctions 0(e) is preserved as well, however, the preservation of 
the realness of fi cannot be concluded based solely on a linear symmetry. 


Assumption (III) (Linear symmetries of A and /). Let A be a densely defined 
and closed operator in a Hilbert space TL and let 5 be a linear , selfadjoint and 
involutive operator, i.e. for all ^ G W and A G C, 5(A0 + 0) = A50 + 50, 
(50,0) = (0,50) and 5 2 = /, such that 
(a) for all 0 G Dom(A), 

50 G Dom(A) and A50 = 5A0, 


(b) for all 0 G A 7 , where y is the subspace from Assumption ! (I) | ([c|), 
5/(0) = /(50) and /(±0) = ±/(0). 


Lemma 4.5. Let A satisfy Assumption EII)1 and let /iq be a simple eigenvalue of 
A. Then 

(i) the eigenvector 0o, ||0o|| = 1, corresponding to the eigenvalue /jlq is either 
S-symmetric or S-antisymmetric, i.e. 

50o = x0o, x = 1 or x = —1. 

(ii) £/ze spectral projection Pq of A corresponding to /io commutes with S, i.e. 


P 0 5 = 5P 0 . 


Proof. The reasoning is analogous to the one in the proof of Lemma 4.2 


□ 


Theorem 4.6. Let A and f satisfy Assumptions ^ T)] an Then the symmetry 

of the eigenvector 0o, corresponding to the simple eigenvalue /io, i.e. 50o = x0o 
with x = 1 or x = — 1, zs preserved for the nonlinear eigenfunction 0 from Theorem 


3.1 i.e. 50 = x0 for all e G (—£q, £q) w tth the same x as for 0q. 


Proof. Similarly to the proof of Theorem |4.4[ we show that 50 = x0 and search 
for the fixed point in a subset B r2£ 2 D {77 G P : Srj = xg}. Note that this is allowed 
since, similarly as in (38)—(39), 


Sx = xx => SG(x’,V = kG(x;<t)- 

5. Applications 


□ 


5.1. Toy model. Let TL = P 2 ((—r,r)) with r = 7r/2, 7 G C and let A 1 be the 
m-sectorial operator associated with the form 

aM := ||V’'I| 2 +7 (IW2)| 2 - |^(-tt/2)| 2 ) , Dom(a 7 ) := H\(- tt/2,tt/2)). 


Since Vi = 0, we take 3^ = tt/2, 7t/2)) = tt/2, tt/2)), c/. Example 2.3 

Note that the inequality, valid for every s > 0, 


< C'llV’IN* ||^>|| < + C(e)||^|| 2 , 


(40) 
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implies that condition ([T]) is satisfied. By standard arguments, cf. |30l Ex. VI.2.16], 
A 7 = —d 2 , 

(41) 

Dom(Ay) = {ip E H 2 ((— 7t/2,7t/ 2)) : 0'(d=7r/2) +7^(±7 t/2) = 0}, 

and (Ay)* = Ay, i.e. A 1 is a non-selfadjoint (unless 7 E M) perturbation of the 
Neumann Laplacian on (— 7 t/ 2, it/2) in boundary conditions; a collection of spectral 
results for A 1 can be found e.g. in [35, [37] . The spectrum of A 1 is discrete with 
explicit eigenvalues: 

<TC) = {-7 2 }U{n 2 }- =1 . 

For 7 = icq a E M \ Z, the operator Ac* is PT-symmetric, cf. [35] and Example 
|4.1[ and all its eigenvalues are real and simple. Hence the spectral condition in 
Assumption (T)](b| is satisfied for any E <7 (Ac*)- The eigenfunctions {£ n }^ 0 5 
{£*}^h 0 °f ^-ia and A* a , respectively, with normalization satisfying Q read 

£o(x) = —j= e~ iax , 

V T 


n > 1, 


[2 n ( / 

( , 7 ry 

\ 1CX • ( 

( , tt\\ 


7 T + - 

-sm 7 

7 T + - 

V 7 T sjn? + a 2 v ' 

V 2) y 

' n V 

V 2)) 


«(*) = 


■£o(A CO) = 


£n(A n > 1. 


sin(<T 7 r) 

We consider the cubic nonlinearity f c from Example |2.8| that is VT -symmetric 
as well, cf. Example |4.1| By Theorems |3.1[ |4.4| and Corollary |3.2[ we have the 
nonlinear eigenpair (/i(e), ip(e)) with ||0(£)|| = 1 and /^(s) €= M for any /iq E <7 (Ac*) 
and small 5. 

For /i 0 = n 2 ,n > 1, straightforward calculations lead to 
t'n = -(Cn|Cn| 2 ,0 = 


and also (p n can be calculated explicitly by solving (22). 

For the special eigenvalue /io = <r 2 and the corresponding eigenfunction 0o = £0, 
the normalized nonlinear eigenpair can be found even explicitly, namely 

M £ ) = a 2 - VY) = ^(0) = 2=e~ iax . 

7r v 77 

Notice the agreement with the expansion of /i and ip from Theorem 


3.1 Indeed, 


since f c (ipo) = A)/ 77 ; identity (21) yields uq = —1/tt. Equation (22) for 0 E 
QoD° m (Ac0 then reads 

—<p>" — a 2 cp = 0, 


hence 0 = 0. Moreover, ip(s) = 0(0) implies that cr = 0 is the fixed point of (31). 


5.2. Bose-Einstein condensates with injection and removal of particles. In 

the physics literature on Bose-Einstein condensates, non-selfadjoint perturbations of 
harmonic oscillators or Laplacians with ^-interactions are considered; the imaginary 
part of the linear potential models the injection and removal of particles, cf. El ED 
l2Tj for instance. 

The nonlinear part of the problem corresponds to the contact (cubic) / c , monopo¬ 
lar / m or dipolar f& interaction, cf. Example |2.10| the unit vector a entering the 
dipolar interaction represents the direction of the polarization, cf. [3HJ EE]. Spec¬ 
tral parameter // corresponds to the chemical potential, the parameter 5 controls 
the strength of the nonlinear interaction and the intensity of particle removal and 
injection (the non-selfadjoint part of the linear operator) is described by the pa¬ 
rameter 7, see (42)—(44). The balance in the removal and injection is reflected in 
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the VT symmetry of the system (the imaginary part of the potential is odd). One¬ 
dimensional examples (obtained by the separation of variables in d = 3 models) of 
“PT-symmetric linear parts A from the literature are the following 


<9 2 + x 2 + v 0 e~ ax2 + i'yxe~ px2 , 

7,^0 Gl,p> 0, 

(42) 

d 2 x + x 2 + vy(<5(x - r) - S(x + r)), 

7 G M, r > 0, 

(43) 

dl ~ (! - h)S{x - t) - (1 + i 7 )< 5 (> + r), 

7 G M, r > 0. 

(44) 


Notice that all can be viewed as holomorphic operator families A(y) with A(0) 
A(0)*, cf. Section 2 . 2.1 


Model (42) corresponds to an operator A from Example 2.2 with cr(A) = crdi SC (A), 
cf. Section 2 . 2.2 For vo and |y| small, all eigenvalues of A are simple and real, cf. Re¬ 
mark 4.3, moreover, it follows e.g. from [42] that the number of non-real eigenvalues 
is finite for any i?o ,7 G M; a numerical analysis of eigenvalues for (42) can be found 
in p T\ . 

Models ([43])—([44]) correspond to the singular Schrodinger operator from Example 
|2.3| Both can be introduced through the closed sectorial form a, cf. ft namely 
Vi(x) = x 2 , v 2 [ip] = i 7 (|' 0 (r )| 2 - |V’(-r)| 2 for ([43]) and Vi{x) = 0, v 2 [ip\ = (-1 + 
i 7 )|' 0 (r )| 2 - (1 + i 7 )|f/>(-r )| 2 fo r (|44| ). Note that condition 0 is satisfied (for any 
7 G C) because of fl40| ). For Q43|), the spectrum is purely discrete and real for 
sufficiently small |y|, cf. Section 2 . 2 . 2 | and Remark |4.3| the number of non-real 
eigenvalues is finite for any 7 G M, cf. [41] for a detailed spectral analysis and m 
for a numerical analysis of eigenvalues. The essential spectrum for (44) is equal to 


[0, + 00 ) and discrete eigenvalues outside of [0, + 00 ) may appear, cf. Section 2 . 2.2 
In more detail, if 


and Remark 4.3 


(1+7 2 )t>1 and (l+y 2 )e 2r > y 2 , 


then there are two negative eigenvalues that are obtained as solutions /i of 

4 


Q —4 


1 


' 7 






cf. [34, 39 . A numerical analysis of eigenvalues of ( [44] ) can be found in [ 8 ]. 

In summary, Theorem |3.1[ Corollary |3.2| and Theorem |4.4| prove the effects ob¬ 
served in physics literature, i.e. for £ 7 ^ 0 , nonlinear eigenvalues are shifted with 
respect to linear ones and those /i that start from real simple linear eigenvalues /io 
are real and the VT symmetry of the nonlinear solution ^ is preserved. Note that 
the latter applies for all nonlinear interactions / po 1 , / m or /d mentioned above and 
a collection of Schrodinger operators in Examples | 2 . 2 | and |2.3| 

A numerical analysis of a model with d = 2 , which is qualitatively similar to 


(42), is performed in Example 6.1 


5.3. Spin-orbit-coupled Bose-Einstein condensate. The spin-orbit-coupled Bose- 
Einstein condensate is described by the spinor <fi G l~i := L 2 (M) x L 2 (M) obeying 
the equation 

id t4> = A(j) + ef((j)) 

with 


f-dl + V{x) + iy uj + iKd x \ 

\ ^ + iftd x — d 2 + V{x) - iyy ’ 


/( 0 ) = 


(( l^l 2 + fe| 2 )M 

V(I^I 2 + I^| 2 W 


(45) 


where V is a trap potential satisfying the conditions in Example |2.2[ k G M is the 
strength of the spin-orbit coupling, w G M is the strength of the linear Zeeman 
coupling, £ characterizes the inter-atomic interactions and 7 > 0 accounts for the 
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decay and the gain of the (pseudo-)spin states up and down, respectively, cf. [29] . 
The time harmonic ansatz </)(x,t) = e~ 1 ^ t 'ip(x) leads to the eigenvalue problem 0 . 

Since the off-diagonal part of A is relatively bounded with respect to the diagonal 
part (i.e. two copies of a Schrodinger operator from Example 2 . 2 ) with the bound 
0, the operator A in ( |45| with Dom(A) = x Hy (M) has the graph norm 

equivalent to W'ipWy 


ll^i 


\n 2 Vi + 11^21 \u 2 v ^ 


X Hy (. 

see Example 


2.2 


Moreover, A has 


compact resolvent if |El(x)| -A oc as \x\ -A oc and it is self-adjoint for real V and 
7 = 0 . Hence, real s imp le eigenvalues (if any) of A for 7 = 0 stay simple and real for 
7 small, cf. Remark [4.3| For instance, if V(x) = x 2 , then, for 7 = 0, the eigenvalues 
of A read 2n + 1 =b uj — a 2 /4, n G No and are all simple if u 7 ^ 0, cf. [29]. Notice 
that the nonlinearity / from ( [45] ) satisfies Assumption |(I)[(|c|) with y := Dom(A), 
|| • ||y as above, cf. the remark on the vector case at the end of Sec. | 2 . 8 | 

If V is PT-symmetric, i.e. V(x) = V(—x), the operator A and the nonlinearity 
/ in (45) possess a natural antilinear symmetry ( C r ip)(x ) := x), where a 1 is 

the Pauli matrix, being the composition of the parity, time and charge symmetries, 
cf. m- Hence our results show the existence of the stationary nonlinear modes 
bifurcating from the linear ones, particularly for the parabolic trap V{x) = x 2 
investigated in |2J|. 


5.4. Optics: nonlinear Schrodinger-type equations. Another set of physical 
applications of nonlinear PT-symmetric problems is optics, cf. for instance IMl 1521, 
where light propagation is modeled by nonlinear Schrodinger-type equations 


i d z u + A u — V (x)u + f(u) = 0 , 


with a gauge invariant /, i.e. f(e ia u ) = e ia f(u ) for all a G M, and typically with 
V G L°°(Q). The variable z is the propagation direction of the optical waves. The 
real part of the optical potential V corresponds to the refractive index and ImV 
models the gain and loss of the medium. If the latter is balanced, in the sense 
that (Imk)(- x) = -(Imk)(r), and if, in addition (ReR)(— x) = (ReH)(x) and 
/ is PT-symmetric, then the whole system is VT- symmetric. For a heterogeneous 
material, V can generally be discontinuous and, in that case, the decomposition to 
Vi and V 2 complying with V\ G must be selected to fit into the setting 

of Example |2. 2 [ however notice that any bounded V fits there after setting V 2 = V. 
The physically most usual nonlinearity is / c , i.e. the cubic one, cf. Example |2.8[ 
Examples of smooth V in d = 2 from m and in d = 1 from jS] are 


V(xi,x 2 ) = — (v 0 + ijX!X 2 )e Xl e * 2 , 

V(x u x 2 ) = -3v 0 ( e -(*i-<d 2 -(* 2 -a ) 2 +e -( Xl +af-(x 2 -a) 2 ' 

— 2 vq ( Xl ~ a ) 2 ~( x 2 -\~ a ) 2 + e -{x!+a) 2 -{x 2 +a ) 2 

— 2 iy ( e -(*i-a) 2 -(* 2 -a ) 2 _ ^—(xi~\-a ) 2 — (x2 —a ) 2 
_ [ry ^ e -(x 1 -a) 2 -(x 2 +a ) 2 _ e -(x 1 +a) 2 -(x 2 +a) 2 ^ 

V(x) = — cos 2 (t) — i7sin(2T), 


7 , v 0 G M, (46) 


7 , t 0 , a G M, (47) 
7 G R. (48) 


2.2.1 


Once again, like for Bose-Einstein condensates, the parameter 7 determines the 
strength of non-selfadjointness, i.e. the loss and gain here, and all models can be 
viewed as holomorphic operator families A( 7 ) with A(0) = 4(0)*, cf. Section 
The potential V satisfies |E(t)| -a 0 as \x\ -A +00 for both ( 
the essential spectrum of corresponding A( 7 ) is [0, 00 ), cf. Section 2 . 2.2 
J R2 ReEdT < 0 (and ReV decays sufficiently fast), there are discrete negative 


and (|47|), hence 
Since 
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eigenvalues of 4(0) = 
tion 


-A + ReK, which are simple for sufficiently small vq , cf. Sec- 


2.2.2 


Neither one of the potentials (46) and (47) is VT -symmetric but both 
are V\T -symmetric and ([46]) is also V 2 T- symmetric, cf. Example 4.1 hence the 


simple real eigenvalues of 4(0) remain simple and real for sufficiently small |y|, 
cf. Remark [4. 3 [ 

Regarding periodic problems, like e.g. ( [48] ) ? our results are relevant for the Bloch 
eigenvalue problem. In the case of 4 = —d 2 + V ( x ) with a 27r-periodic V, one 
considers the family of operators 4 in L 2 ((— 7 r, 7 r)) with /c-quasi-periodic boundary 
condition, i.e. with the domain 

Dom(4) = {'ip E H 2 ((— 7 r, 7 r)) : 0( tt ) = e l/c 0(— 7 r), ^'(tt) = e lk/ ip'(— 7 r)} 


2.5 


Since V in (1481) is VT- 


and the form domain l~L]f k ((— 7 r, 7 r)), cf. Example 
symmetric, eigenvalues of A for k 0 {0, tt} are simple and real for sufficiently small 
|y|, cf. Section 2 . 2.2 and Remark [43} Numerical analysis from [44] suggests that 


for (48) this is the case if |y| < 1/2. 

In summary, our results in Theorem |3.1 
applicable and provide for 


Corollary |3.2| and Theorem |4.4| are 
( |46| ) and (47) and any / compatible with Assumptions 
|(I) | and [(II) | real nonlinear eigenvalues p of — Aip + Vip — sf(^p) = /i 0 with nonlinear 
solutions 0 that satisfy the corresponding partial VT symmetries. For the periodic 
problem ( |48| ), we obtain nonlinear Bloch functions ip(x) = p(x)e lkx , where p is 2tt- 
periodic. This complements the results of [18] on the bifurcation of nonlinear Bloch 
waves in the selfadjoint case. In the ^-dependent nonlinear Schrodinger equation 
we obtain solutions u(z,x) = e~ l ^ z ^(x) with a real propagation constant p despite 
the fact that the material exhibits loss and gain. 

A numerical analysis of the model with the potential in (47) is performed in 
Example | 6 . 2 | 


5.5. Optics: discrete nonlinear Schrodinger equation. The propagation of 
light in a finite one dimensional lattice of linearly coupled Kerr-nonlinear fibers is 
often modeled by the discrete nonlinear Schrodinger equation 


id z u n — ^n+l T 1 T h/C 1) ^ n T \u n 


1 < n < 2 A, uq = R 2 A +1 = 0 , 


where z is the propagation direction, n E N is the lattice site and ry(—l) n E iM 
describes the loss or gain at the site n, see m- For time harmonic solutions 
u n (t) = e _ 1 Att 0 n and after the rescaling ip n ::= e _ 1 / 2 0 n (with e > 0 ), we get eigen¬ 
value problem 0 with A in ( [To] ) and the nonlinearity / n ( 0 ) = IV’nPV’n- Exam¬ 
ple [TGI and Section |2.2.3 explain that for |y| small enough Assumption m holds 
oj — m — rp 2 ? 7 V Note that the Lipschitz continuity of / holds, e.g. with 


with T-L = y 

IV 


2^7 = 1 


C 2 

|071 


2 0n | < CmaXj = i j ... j2 iv{|'0n 


,|0n| 4 }E-fll0n 
(i.e. the discrete PT-symmetry) 
such that V- n = V n . Our 


Assumption (II) is satisfied with (C0) n = 
due to the choice of the “potential” V n := ry(—1) 
results therefore recover Theorem 1 in [31]. 


5.6. Optics: coupled mode equations. In Kerr-nonlinear optical fibers with 
a Bragg grating and a localized defect the propagation of asymptotically broad 
wavepackets can be described by the system of “coupled mode equations” 

i (dtEh. + d x E{) + k{x)E 2 + V(x)E 2 + (l^l 2 + 2|£ 2 | 2 )£ 1 = 0 
i (d t E 2 - d x E 2 ) + k(x)Ei + V{x)E 1 + (|£ 2 | 2 + 2\E 1 \ 2 )E 2 = 0 

with k(x) -a Kqq > 0 and V(x) -A 0 as \x\ -A oc, see [22]. The potentials k(x) — Koo 
and V ( x ) describe the defect of the material and are determined by the refractive 
index. Once again, we consider the time harmonic ansatz E(x,t) = e~ 1 ^ t (p(x) and 
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after the rescaling i/j := e x ! 2 (\) (with 6 > 0), we obtain eigenvalue problem 0 with 


A in (11) and 


fh/j) -f( IV’i| 2 + 2|V’ 2 | 2 )A 
fW ~ U | V ’ 2| 2 + 2 | V > 1 | 2 )^2 


As mentioned in Section [2. 2. 4[ real smooth and bounded potentials a and V exist 
such that A has a simple isolated eigenvalue. Examples |2.7| and |2.8| (see the remark 
on the vector case at the end of Sec. 2.8), guarantee that Assumption | (I) | is satisfied 
with U = L 2 (R) x L 2 (R) and y = H\R) x H\R ) provided V, k G L°° and 
k(x) -a Koo > 0 and V(x) -A 0 as \x\ -A oo. 

For materials with loss/gain the potentials V and a become complex and choos¬ 
ing them PT-symmetric, we satisfy also Assumption |(II)| The existence of a real 
simple isolated eigenvalue /io of A is guaranteed at least for small imaginary parts 
of k and V by the analytic dependence as in Remark |4.3| Hence, (in the language 
of [22]). our results show that conservative nonlinear defect modes bifurcate from 
linear ones in the VT -symmetric case. 


5.7. Superconductivity. A model of a finite superconducting wire is discussed in 
[48l 49 and the bifurcation of nonlinear states for a nonlinear parabolic equation 
(d = 1) on Q = (—1,1) is studied. In detail, the problem 

w t = w xx + i xlw + Tw + N[w\) x G (—1,1), w(—l) = w( 1) = 0, 

f x (49) 

N[w] = — \w\ 2 w + iw / Im (w(s, t)w x (s, t)) ds, 

Jo 

where / and T are real parameters, is considered. In [49] Sec.6] the authors study 
the bifurcation of nonlinear (generally t-dependent) solutions from the zero solution 
at the smallest eigenvalue Ai G R of 

A:=-d 2 x -ixI, Dom(A) :=ff 2 ((-l,l))nff 0 1 ((-l,l)). (50) 

The potential —i xl is PT-symmetric, so the spectrum of A remains real if the 
parameter I is chosen small enough, cf. Remark |4.3[ and the number of non-real 
eigenvalues remains finite for any I G R. For the bifurcation problem the authors 
set T = ReAi +£,0<£<1 and use the center manifold reduction, where the center 
manifold is one dimensional and corresponds to the zero eigenvalue of A — ReAi. 
On the manifold they study ^-dependent, but also stationary nonlinear solutions. 
The asymptotics of the latter are given by 

w{pc) ~ e 1 ' 2 aui(x), 

where u% is the linear eigenfunction corresponding to Ai and a G R is the pro¬ 
jection coefficient on the center subspace and solves an algebraic equation. For 
I small enough Ai G M, such that a real nonlinear eigenvalue T bifurcates. The 
eigenfunction w is VT -symmetric due to the VT -invariance of the center manifold. 
In the formal part of m the more detailed expansion 

w(x) ~ s 1 / 2 aui(x) + s 3 / 2 wi(x), 

is given, where the correction w\ solves 

(A — \i)wi = aui + N\aui]. (51) 

a G R can then be selected via the solvability condition of the above equation and 
agrees to leading order with the a from the center manifold approach. 

To relate this work to our results, we rescale the ^-independent solution w(x) = 
e 1 ^ 2 2 p(x) and recover from ( |49| ) a problem of type 0, namely 

{A — r)t/> — e(/ c ( , 0) + /n(V’)) = 0 , 
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cf. Examples 2.8, 2.11 Equation (51) thus corresponds to our ( 22 ). Observe 
that A fits into the setting of Example 2.3 with V\(x) = ix/, V 2 = 0 and y = 
Hn’ Dir ((—1,1)) = Hq((— 1,1)). The nonlinearities are discussed in Examples 


2.11 


and |4.l| and it shown that H 1 is a suitable space for the Lipschitz condition (|4|. 


Hence, our results in Theorem |3.1| Corollary |3.2| and Theorem |4.4| are applicable 
and provide real nonlinear eigenvalues /x with PP-symmetric nonlinear solutions p 


6. Numerical Examples 


We analyze numerically two nonlinear problems of type (HI), both with the 
Schrodinger operator A = —A + V in L 2 ( 1 R 2 ), cf. Example |T 2 [ and the cubic 
nonlinearity / c , cf. Example |2.8[ i.e. 

(~A + V)ip - e\ip\ 2 ip = mV’, IIV’IIl 2 = 1- (52) 


Selected potentials V posses some antilinear or linear symmetries. Clearly, the 
nonlinearity is highly symmetric and satisfies Assumptio n |(II) (|b|) with C = VT 
as well as C — VjT , j = 1 , 2 , and also Assumption (III) (|b| with any coordinate 
reflection symmetry <S. 

Our choice of d = 2 rather than the numerically simpler d = 1 allows the 
investigation of partial PT-symmetries as well as the interplay between antlinear 
and linear symmetries in a single problem. 

The numerics are performed using the package pde2path Hannon for numeri¬ 
cal continuation and bifurcation in nonlinear elliptic systems of PDEs. The package 
uses linear finite elements for the discretization, Newton’s iteration for the compu¬ 
tation of nonlinear solutions and an arclength continuation of solution branches. 
In all numerical computations, the free complex phase of the solution was fixed by 
forcing Im^O^o)) = 0 at a selected point xq within the computational domain. 
For the plots, we select xq = (0,0) for the PT-symmetric case in Example 6.1 and 
xq = ( 0 , 2 ) for the PiP-symmetric Example 6.2 In all computations, except for 


one case mentioned below, the numerical grid is selected symmetric with respect 
to both coordinate axes as well as with respect to the reflection x —x. This is 
crucial for recovering symmetries of eigenfunctions and realness of eigenvalues. 


Example 6.1. We consider first the following imaginary perturbation of the har¬ 
monic oscillator that is compatible with Example | 2 . 2 | and inspired by the Bose- 
Einstein condensates models from Section [572] , 

V(x 1 ,x 2 ) = + x t) + ^2 + ^2 + 2 ’ ( 53 ) 

Clearly, V(—xi, —X 2 ) = V{x\^X 2 ) = V(— aq,^)- Hence, the problem has three 
symmetries: two antilinear symmetries, namely the VT symmetry and the Pi T 
symmetry, and the linear V 2 symmetry, cf. Example |4.1| 

For 7 = 0, the eigenvalues of A are known explicitly: 

A/e = \/2k, k = 1 , 2 ,..., where the multiplicity of A/e is k. 

Enumerating the eigenvalues including their multiplicity, we obtain our eigenvalues 
H n for 6 = 7 = 0 . 

For the discretization of the PDE, we take 2*80 2 = 12800 isosceles right triangles 
of equal size generated by Matlab’s command poimesh on the domain x G [— 8 , 8] 2 
with homogeneous Dirichlet boundary conditions. The first four eigenfunctions are 
well localized within the selected domain. 

For 7 = 2, the numerically obtained first four eigenvalues (for 6 = 0) are 

fjbi ss 2.096, /i 2 ~ 2.583, /x 3 ss 3.155, /X 4 ~ 4.256, 
and they are all simple. 
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(a) 7=0 (b)7=2 




Figure 1. Bifurcation diagram of (|52j), (|53 
for the first eigenvalues /ii, 


in the parameter s 
/i 4 with 7 = 0 in (a) and with 7 = 2 


in (b). Circles label secondary bifurcation points. 


In Fig. |6.1| the bifurcation diagram in 5 for fi i,/i 2 ,M 3 and /i 4 is plotted. The 
eigenvalues are continued in 6 > 0 from the linear eigenvalues at e = 0 for two values 
of 7 , namely 7 = 0 in (a) and 7 = 2 in (b). In both cases, all plotted eigenvalues 
(including the ones bifurcating from ^ 3 ) are real. Note that for 7 = 0 the grid is 
symmetric only with respect to the reflection x -A —x and not with respect to the 
coordinate axes. This suppresses the multiplicity of the first four linear eigenvalues 
such that these can be easily numerically continued in 6 . The numerics suggest 
that all the four simple eigenvalues remain real for at least e < 10. Clearly, the 
numerics agree with the analysis as simple eigenvalues stay real for e small. In 
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Figure 2 . Profiles of the nonlinear eigenfunctions of (52), (53) at 
7 = 0 labeled by 1-4, 3b in Fig. |6.1| (a). 


Fig. [2J the eigenfunctions at the five labeled points at e = 5 in Fig. | 6 . 1 | (a) for 
7 = 0 are plotted. Since V is real, the eigenfunctions can be automatically chosen 
“PT-symmetric (using a proper rotation of the complex phase). Note also that after 
a proper rotation in the (xi, X 2 )— plane (allowed due to the rotation symmetry of 
V at 7 = 0 ) all eigenfunctions 1-4 and 3b are symmetric or antisymmetric with 
respect to x\ —x\ as well as X 2 —X 2 - The numerically generated profiles for 
the eigenfunctions 1-4 are symmetric about other axes due to the lack of coordinate 
symmetry of the grid, as explained above. 

In Fig. [3J the profiles for the case 7 = 2 from Fig. |6.1| (b) appear. Eigenfunc¬ 
tions 1 and 2 satisfy all the three symmetries, i.e. VT, V\T as well as the linear 
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Figure 3. Profiles of the nonlinear eigenfunctions of (52), (53) at 
7 = 2 labeled by 1-4, 3b in Fig. |6.1| (b). 


^(aq, —X 2 ) = ±^( 2 ; 1 ,^ 2 )- Eigenfunctions 3 and 4 satisfy the linear (anti)symmetry 
and can be chosen either VT- or V\T -symmetric after a suitable multiplication by 
i. The eigenfunction 3b on the dotted branch (bifurcating from the primary branch) 
is only V\T- symmetric. 

In Fig. [4j we perform continuation in the parameter 7 for the two values 5 = 0 
and 5 = 2. The results are qualitatively similar to those in ID from El. When two 
real eigenvalues collide, they leave the real axis and become a complex conjugate 
pair. In addition, however, a secondary bifurcation can occur, like, e.g., from fi 3 at 
5 ss 1.5 (see the inset in Fig. 0 (c)). 

Note that for a complex conjugate pair of simple eigenvalues the correspond¬ 
ing eigenfunctions ^ 1^2 can be chosen to be related by ^{x) = (VT'ipi)(x) = 
ol(ViT^i){x) with a suitable factor a G C, \a\ = 1. Hence, below we always plot 
only one eigenfunction for a complex conjugate pair. 

The eigenfunctions for the linear case 5 = 0 at the six labeled points in Fig.^a) 
are shown in Fig. [5j The symmetry properties of the eigenfunctions 1-4 are the same 
as for the case 7 = 2,5 = 5 in Fig. [ 3 ] The eigenfunctions 5 and 6 corresponding 
to the complex eigenvalues are neither VT- nor TiT-symmetric but the linear 
symmetry ^(^ 1 ,-^ 2 ) = ±^(^ 1 , ^ 2 ) is preserved. This is in agreement with Lemmas 
l4~2l and I3~5l 

For the nonlinear case 5 = 2, the eigenfunctions are in Fig. [ 6 j The symmetries 
of the eigenfunctions 1-4 are again the same as for the case 7 = 2 ,5 = 5 in Fig. [3j 
The eigenfunctions 5 ad 6 of the complex eigenvalues are, once again, only linearly 
symmetric. On the other hand, the eigenfunction 3b on the secondary bifurcation 
branch (corresponding to a real eigenvalue) has only the V\T symmetry. 


Example 6.2. As the second example, we choose the T^iT-symmetric potential 
© with a = | and v$ = 1 , appearing in optics literature. Clearly, V satisfies 
V(— ^ 1 ,^ 2 ) = F(xi,£ 2 ), but no other obvious antilinear or linear symmetry involv¬ 
ing reflections of coordinates. Simplicity of eigenvalues of A is discussed in Section 


5.4; the four lowest eigenvalues appear to be simple numerically for 7 = 0 . 


Our discretization mesh is given by 14400 isosceles right triangles on the domain 
x G [—13,13] 2 chosen such that the mesh is symmetric about the coordinate axes 
and with respect to the reflection x -A — x. We again use homogeneous Dirichlet 
boundary conditions. 

In Fig. [7] we plot the bifurcation diagram in the parameter 5 for 7 = 0 and 
7 = 0.1. Similarly to Example 6.1 for 7 = 0 (real potential F), the nonlinear 
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Figure 4. Bifurcation diagram of (52), (53) in the parameter 7 
for the first eigenvalues /xi,... , 7 x 4 with £ = 0 in (a) and (b) and 
with e = 2 in (c) and (d). Circles label secondary bifurcation 
points. 
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Figure 5. Profiles of the nonlinear eigenfunctions of (52), (53) at 
£ = 0 ,7 = 1 labeled by 1-6 in Fig. [4] (a). 


eigenvalues stay real even after secondary bifurcations, while for 7 = 0.1 secondary 
bifurcations result in complex conjugate pairs of nonlinear eigenvalues. 

In Fig.[ 8 j we plot the six eigenfunctions labeled in Fig. 0 (b) for the case 7 = 0 . 1 . 
Clearly, all the four eigenfunctions on the primary branches (labels 1-4) with real 
eigenvalues are PiT-symmetric. The bifurcating solutions (labels lb and 3b) are 
asymmetric. Once again, for complex conjugate pairs we plot only one eigenfunction 
as the two can be chosen to be related by ^2 (x) = 

The bifurcation diagram in the parameter 7 is plotted in Fig. [9] for the two 
values £ = 0 and £ = 2 and the results are, again, analogous to Example |6.1| The 
first collision of eigenvalues occurs at 7 = 7 * ~ 0.22 in agreement with the value 
7 * = 0.214 reported in m • 
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Figure 6 . Profiles of the nonlinear eigenfunctions of (52), (53) at 
5 = 2 ,7 = 1 labeled by 1-6, 3b in Fig. [4] (c). 


Figures [lO] and [TT] show the eigenfunctions at the points labeled in the bifurcation 
diagram in Fig. [9] The eigenfunctions at points 1-4 (before eigenvalue collision) are 
all V\T -symmetric while after the collision at points lb and 3b they are asymmetric. 
This is in contrast with Example |6.1[ where linear symmetry was preserved in 
collisions. Here, no obvious linear symmetry is available. 
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Figure 7. Bifurcation diagram in the parameter 6 for the first 
eigenvalues ft i,... ,/m of ( [52] ) ? ( |47| ) with 7 = 0,uo = 1 in (a) and 
with 7 = 0.1, no = 1 in (b) and (c). Circles label secondary bifur¬ 
cation points. 
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Figure 9. Bifurcation diagram in the parameter 7 for the first 
eigenvalues /xi,..., /14 of of ( [52] ) , ( |47| ) with 5 = 0, vq = 1 in (a) and 
(b) and with 6 = 2,uo = 1 in (c) and (d). Circles label secondary 
bifurcation points. 
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